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We perform a comparative theoretical study of the data at spacelike momentum transfer for the 
7*7 — > 7T° transition form factor, just reported by the Belle Collaboration, vs. those published before 
by BaBar, also including the older CLEO and CELLO data. Various implications for the structure 
of the 7T° distribution amplitude vis-a-vis those data are discussed and the existing theoretical 
predictions are classified into three distinct categories. We argue that the actual bifurcation of the 
data with antithetic trends is artificial and reason that the Belle data are the better option. 
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The newly released data for the pion-photon tran- 



sition form factor (TFF) F~« 77r (Q 2 ,q 2 -> 0) by the 
Belle Collaboration [l| , seem to be dramatically different 
from those reported in 2009 by the BABAR Collaboration 
Instead of a pronounced growth of the scaled TFF 
with Q 2 , observed by BABAR above 9 GcV 2 , the Belle 
data are compatible (with the exception of one point) 
at high Q 2 with the limit set by perturbative QCD Q: 
q2 F i'W = « 0.185 GeV with f n = 131 McV. 

The behavior of the tt° TFF with increasing Q 2 can 
be forecast within calculable theoretical errors using the 
theoretical tools of QCD — perturbative and nonpertur- 
bative. Perturbative QCD governs evolution and supplies 
the means to calculate radiative corrections, whereas 
nonperturbative QCD provides the modeling concepts 
and techniques to determine the ir° distribution ampli- 
tude (DA) — in leading-twist two and higher twists — and 
treat the hadronic structure of the nearly on-mass-shcll 
photon. 

In a recent paper we have discussed what one 
should expect for the outcome of a measurement of 
the 7T° TFF at high Q 2 according to QCD. We argued 
that the expected scaling behavior should follow — within 
uncertainties — the interpolation formula of Brodsky- 
Lepage g: F^(Q 2 ) = (y/2f„) / (4tt 2 / 2 + Q 2 ). This 
phenomenological formula links the TFF at Q 2 = 0, fixed 
by the axial anomaly, with the QCD asymptotic limit 
\2fw While plausible and useful in practical terms, 
this formula is not derived from QCD. Hence, it is of 
paramount importance to calculate the TFF within QCD 
in order to obtain an expression that can replicate the 
appearance of stasis in the scaled TFF above some Q 2 
value as a result of switching on parton-photon interac- 
tions controlled by QCD. Such a saturating behavior of 
the TFF can be obtained from a formalism developed 
in a series of papers [rljioj that combines the dispersive 
approach of light-cone sum rules (LCSR)s [III 12 1 ( see 
also 13, |l4j) with QCD sum rules that employ nonlocal 
condensates [ljj [l6[. The latter scheme is used to de- 
rive the pion DA, while the former serves to properly 



accommodate the hadronic content of the low-virtuality 
photon. Then, the TFF is defined by the LCSR 
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given by 

pert(n2 



where the spectral density 
p(Q 2 ,x) = (Q 2 + s)pP crt (Q 2 ,s) with pP crt (Q 2 ,s) 
(l/7r)ImF 7 7 * (Q 2 , — s — ie), and the abbreviations 
x = 1 — x, s = xQ 2 /x, xo = Q 2 /(Q 2 + s ) have been 
used. The hadronic content of the quasi-real photon 
in the TFF is taken care of by means of the first term 
in Eq. (fTJ), whereas the partonic pointlike interactions 
above s > sq are described by the second term which is 
calculable order by order in QCD perturbation theory. 
Our computation is performed at the level of the NLO 
spectral density by taking into account the correction 
pointed out i n Il4| . The other parameters have the fol- 
lowing values 0: s = 1.5 GeV 2 , m p = 0.77 GcV, while 



the Borel parameter is M 



Ml pt /(x) 



< 1 GcV 2 



from the two-point QCD sum rule for the p-meson with 
G [0.5 -7- 0.8] GeV 2 (see for more details). 
The twist-two pion DA in the formalism of Bakulev, 
Mikhailov, and Stefanis (BMS) @ is modeled in terms 
of two Gegenbauer coefficients a 2 and 04: ip^ BMS (x) = 



^ sy {x) 
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where 



</? asy (a;) = 6xx denotes the asymptotic (asy) 7r° DA. 
The values of a 2 (^ 2 ) = 0.20 and a 4 (^ 2 ) = -0.14 (at 
/i 2 = 1 GeV 2 ) are selected in such a way that the first 
ten moments (^ 7V ) 7r = Jq dx(2x — l) N tp^k\x, fi 2 ) with 
the normalization condition dxip^\x, u 2 ) = 1 lie 
inside a particular range determined in [6|, while all 
higher coefficients ae,as,aio were determined and found 
to be negligible. This procedure gives rise to a whole 
"bunch" of possible DAs whose shape is characterized 
by a double-humped structure with strongly suppressed 
endpoints x = 0, 1. This suppression is controlled by the 



2 



vacuum quark virtuality X 2 q (fi 2 = 1 GcV 2 ) « 0.4 GcV 2 
which characterizes the nonlocality of the quark (quark- 
gluon) condensate [f| and corresponds to a correlation 
length of about 0.31 fm. This parameter turns out to 
be intimately related to the Twist- four coupling S 2 (fj, 2 ) 
which has a value around S 2 ps X 2 / 2 — details can be 
found in 0. The BABAR. data from Q 2 = 10 GeV 2 
onward show a marked tendency to increase with Q 2 
and are therefore incompatible with the BMS formalism. 
These data can be best reproduced with a flat-top 7r° 
DA [li-iOl that yields an auxctic TFF 0— more below. 

Let us look more closely in Fig. Q] at what the exper- 
imental data from different collaborations CELLO |2l| . 
CLEO [13], BABAR 0, and Belle [l| mean, comparing 
them with the results of several theoretical approaches. 
The shaded (green) band shows the predictions calcu- 
lated using the BMS formalism detailed in [l(|. The re- 
sult for the BMS model @ is represented by the solid line 
inside it, while the width of the band collects uncertain- 
ties from different sources: (i) the variation of the shape 
of the 7T° DA extracted from QCD sum rules with nonlo- 
cal condensates @ , (ii) the uncertainty of the Twist-four 
coupling S 2 0, and (iii) the sum of the Twist-six term 
[lH and the next-to-next-to-leading order (NNLO) ra- 
diative correction, proportional to the /3q function, com- 
puted in |22j . The combined treatment of these last two 
uncertainties is justified because for the Borel parameter 
M 2 < 1 GeV 2 , we are adopting, see, e.g., [ijj], both have 
rather small magnitudes comparable in size but opposite 
signs, with the Twist-six term being positive. 

Dropping extraneous details, let us address the other 
curves shown in this figure from bottom to top. The 
dashed line at the lower border of the shaded (green) 
band — tagged Asy — denotes the TFF computed with 
the asymptotic n° DA. Next, the dotted [23| and the 
double-dotted-dashed [24| lines, partly intersecting with 
the band, give the results obtained with two holographic 
models based on the AdS/CFT correspondence. The 
thick dashed-dotted line at the upper boundary of the 
shaded band shows the recent prediction for the TFF ob- 
tained in [25} using an extended vector-dominance model. 
The intermediary (blue) solid lines denote the TFF ob- 
tained for models III (lower curve) and I (higher curve) 
from the LCSR analysis in [13], whereas the prediction of 
their model II lies in between (not shown) . The key char- 
acteristic of these models for the 7r° DA is the large size 
of the coefficient 04 > 02 and the non-negligible values of 
the higher Gegenbauer coefficients a^, a 8 , a\o, ai2- Actu- 
ally, model I is a flat-top DA with a reduced coefficient 

a flat-to P( y = X Gc y2) _^ q 130) model m 

coefficients 02 = 0.160, 04 = 0.220, a 6 = 0.080 (all values 
at 1 GeV 2 ) [lH — neglecting higher ones. Note that this 
analysis takes into account the twist-six contribution — 
found to be positive and very small for the adopted value 
of the Borel parameter M 2 = 1.5 GeV 2 , while it ignores 
the explicit inclusion of the (negative) NNLO radiative 



correction. As one sees from this figure, the inverse hier- 
archy of the Geg enbauer coefficients a-i < 04 in the LCSR 
approach of [14j turns out to be not really enough to fully 
account for the auxetic behavior of the BABAR data above 
9 GeV 2 , while for the same reason it causes a deviation 
from the Belle data and the asymptotic QCD limit. 

A similar result — shown as a double-dotted-dashed 
(blue) line (Fig. [TJ) in the vicinity of the previous two 
curves — was obtained by Kroll in [26| using a different 
theoretical framework based on the Botts-Li-Sterman 
Modified Factorization Scheme 27, Eil . At the upper 



end of this regime of predictions one finds the TFF (long- 
dashed-dotted line), computed by Polyakov in |20| by 
using a flat-like tt° DA, that was extracted from an effec- 
tive chiral quark model based on the instanton vacuum 
of QCD and including leading-order (LO) evolution ef- 
fects. The two (red) curves shown at the top represent 
the predictions obtained with a constant flat-top DA by 
Radyushkin [lj| (dashed line) and Lih and Geng [29[ 
(thick dashed-dotted line). For completion, also the pre- 
diction for the TFF for the Chernyak-Zhitnitsky (CZ) tt° 
DA [3(| is shown (solitary line with the CZ label). 




FIG. 1: (color online). Theoretical predictions for the scaled 
7*77T° TFF, calculated in various approaches, in comparison 
with all existing experimental data, using for the latter the 
indicated notations. The shaded (green) band contains the 
results obtained within the BMS formalism in [l(j, with the 
solid line inside it denoting the prediction for the BMS model 
DA The other curves are explained in the text. 

The main message from Fig. Q] is that only with a flat- 
top pion DA one can really replicate the auxetic behavior 
of the 7T° — > 7*7 transition form factor exhibited by the 
high-Q 2 BABAR data [l8l420l |. It is worth noting, however, 
that to get best agreement in the statistical sense, one has 
to "cull" the two outliers at 14 and 27 GeV 2 before analy- 
sis. These two outliers are close to the asymptotic QCD 
limit and hence, strictly speaking, incompatible with a 
flat-top 7T° DA. From the theoretical point of view, a flat- 
top DA reflects a contingency approach geared in hind- 
sight with the only aim to emulate the auxetic trend of 
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FIG. 2: (color online) . Classification of theoretical predictions 
for the spacelike scaled TFFs of the tv° and the nonstrange 
component |n) of the 77, and 7/ in terms of three distinct 
shaded bands in comparison with all existing experimental 
data with notations as indicated. The predictions for the 
asymptotic (lowest dashed line) and the CZ DAs (solid line) 
are explicitly labeled. 



the BABAR data for the pion TFF. There are no features 
in some region of the longitudinal momentum fraction x 
of the valence quark and no humps anywhere. Of course, 
this scaling must eventually break down at the endpoints 
x = and x = 1 because no parton can carry a larger 
momentum fraction than 1. Therefore, a flat-top ir° DA 
describes the pion as being a "pointlike" particle with no 
internal structure because it looks the same everywhere 
j3ll |. On the other hand, the BMS formalism yields pre- 
dictions that are possible only if the ir° TFF is governed 
by QCD — in lingo: collincar factorization underlying Eq. 
([1]). Accordingly, the TFF saturates at high momentum 
values of Q 2 , meaning that the high-virtuality photon 
couples to a single parton inside the considered meson 
which is well-described by a BMS-like DA that incorpo- 
rates the nonlocality of the QCD vacuum. In contrast, 
the BABAR data Q for the pion contradict this behavior 
necessitating some nonperturbative mechanism capable 
of providing, e.g., ln(<5 2 //x 2 ) enhancement (where [i 2 is 
some contextual nonperturbative scale < 1 GeV 2 ) to the 
TFF, as encoded in a fiat-top tt° DA (H-E^. In short, 
the invention of a flat-like DA within a specific context is 
contingent on the unforeseen behavior of the n° BABAR 
data at high Q 2 . Would the outcome of these data be 
in line with the asymptotic QCD prediction, attempts to 
create a flat-like w° DA to challenge the collinear QCD 
factorization would appear contrived and artificial. Re- 
lying instead solely on the Belle data [l| , there would be 
no need at all to invoke a fiat-like DA based explanation 
of the 7T° TFF. 

The dichotomy between scaling of the TFF with Q 2 
and auxesis becomes more evident by inspecting Fig. [2] 
which summarizes the experimental situation for the tt° 
TFF in conjunction with the predictions of various theo- 



retical approaches. The latter can be organized into three 
distinct shaded bands, whose widths are adjusted to in- 
clude similar predictions obtained in different theoretical 
schemes and are not related to particular inherent the- 
oretical uncertainties (like the BMS band shown in Fig. 
[1} . The main goal of presenting Fig. [2] is to place the 
7T° TFF in a broader perspective by collecting and com- 
paring what is known experimentally and categorizing 
what has been proposed theoretically. The central obser- 
vation is that the experimental data for the ir° TFF ar- 
range themselves at high-Q 2 values (starting at 10 GeV 2 ) 
along two branches which below this value unify into a 
single one. The upper branch of the data consists only 
of the BABAR data for the tt° TFF 0]— with the excep- 
tion of the two outliers downwards at 14 and 27 GeV 2 , 
already mentioned in the discussion of Fig. [TJ and the 
Belle outlier upwards at Q 2 = 27.33 GeV 2 . The lower 
branch contains the two BABAR outliers and all the other 
Belle data. One can also count to this branch the data 



ex- 



4] 



diere \n) = (1/2) (\uu) + \dd)), 



for (3/5)Q 2 Ft 1 

tracted from the BABAR data (32| for r\ and rf — see 
for more. 

The classification of the theoretical predictions fol- 
lows roughly this pattern plus an additional band in 
between. Those predictions agreeing with the standard 
QCD factorization scheme are forming the lower (dark- 
green) band, labeled 1. To be specific, this band contains 
the predictions obtained within the BMS formalism in 
our most recent analysis in fiol 33-35J, where one can 
find the details. Moreover, band 1 includes the result of 
the form- factor modeling of (25[, based on an extended 
vector-dominance model. Interestingly, also the predic- 
tions from two different AdS/QCD models, viz., [24[ and 



23l l36l |. lie within band 1. The lower boundary of this 
band is compatible with the asymptotic tt° DA (dashed 
line with the flag Asy). 

The (red) band 3 collects the results from [lj| and 
j29l ] , which both employ a flat-top 7r° DA, as well as that 
of the analysis in [37j which utilizes a flat-like DA and 

3£ 



and \4£ 



Sudakov effects. The predictions from [38 . 
based on a dispersive representation of the axial anomaly 
and quark-hadron duality, are within this band as well. 
The similar results of jUJ are also incorporated, while the 
findings of 42, 43| (not shown) would appear just below 
the edge of band 3. Moreover, band 3 contains the result 
of the calculation in [44| which ascribes the auxetic TFF 
behavior to new physics in the r sector. 

Between the two aforementioned bands, one has an- 
other class of theoretical predictions forming the (blue) 
band 2. This consists of the results obtained with mod- 
els I, II, and III of the LCSR analysis of [l4[ , Polyakov's 
result [20j, extracted from the chiral quark model, and 
also an analogous result from [jjj . Note that the curve 
representing the TFF computed with the CZ DA (shown 
as a single solid line with the label CZ) joins band 2 at 
the far end of the experimentally accessible momentum 
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TABLE I: We show in the first row the xV n df for the BMS model in comparison with estimates of the coefficients a n of 
the 7T DA determined by fitting the pion-photon TFF within LCSRs. The second and third rows show a 2D fit in the (02, 04) 
plane, while the last two rows employ a nonzero coefficient a^,. The errors are due to statistical uncertainties and a systematic 
error related to the Twist-four term. The last column shows the values of "x 2 /mil (with ndf = number of degrees of freedom) 
for the considered data sets. [All entries evaluated at /i| Y = 5.76 GeV 2 with SY abbreviating Schmedding and Yakovlev 



Data set 


«2(msy) 


04(Msy) 




X 2 /ndf 


CELLO [21], CLEO [17], Belle [1] 


0.142 


-0.090 





22.1/33 


CELLO, CLEO, Belle 


0.154 ± 0.046 ±0.055 


-0.066 ± 0.067 ± 0.058 





20.1/31 


CELLO, CLEO, RABAR [2] 


0.090 ± 0.037 ±0.050 


0.069 ± 0.057 ± 0.053 





69.5/33 


CELLO, CLEO, Belle 


0.157 ±0.057 ±0.056 


-0.192 ±0.122 ±0.077 


0.226 ±0.161 ±0.033 


13.1/30 


CELLO, CLEO, BABAR 


0.177 ± 0.054 ±0.056 


-0.171 ±0.103 ±0.071 


0.307 ±0.096 ±0.024 


33.3/32 



region, while it is more than 4tr away from all data below 
10 GeV 2 . Except the calculations giving rise to strip 3, 
most other predictions shown in Fig. [5] include Efremov- 
Radyushkin-Brodsky-Lepage (ERBL) 0, 47 1 evolution 
at the LO or NLO level. 

In Table HI we present statistical fits of the pion TFF, 
calculated with LCSRs to the data. The first row shows 
the original BMS values from [(| , while the next two rows 
show a 2D fit to two sets of data using only the first two 
coefficients 02 and 04. One set contains the CELLO, 
CLEO and Belle data, while the other one consists of 
the CELLO, CLEO and BaBAR data. The last two rows 
represent an analogous fit which also employs the next 
coefficient ae. Evidently, the Belle data allow in both 
cases a better statistical description. This distinct be- 
havior completely agrees with our classification scheme 
presented in Fig. [21 with the BMS model and the 2D 
best-fit being entirely inside band 1. 




FIG. 3: (color online). Unified range of theoretical predic- 
tions (grey area) discussed in context with Fig. [2] The darker 
shading indicates transition from scaling with Q 2 to auxesis. 
The experimental data of CLEO, Belle, and BABAR are shown 
in terms of fits (single lines) explained in the text. 

But from another more data-oriented point of view, 
one may argue (see, e.g., 



tive deviation between the Belle and the BABAR fits does 
not exceed 1. 5a [l[. Appealing to this comparability, one 
may be tempted to unify the theoretical results, shown in 
terms of three distinct strips in the previous figure, into 
a single wide band. This is illustrated in Fig. [3J where 
the darker shading of the band towards the top indicates 
the deviation from scaling with Q 2 predicted by QCD 
towards an auxetic TFF behavior contingent on model- 
dependent explanations. Thus, even if such data pooling 
seems statistically acceptable, the underlying theoretical 
approaches are hardly comparable to each other, so that 
an interpretation of the Belle and the BABAR data in re- 
lation to each other and against some common standard 
appears rather questionable. In this figure, all data sets 
are represented by single fits as follows: (i) The top solid 
line shows a power-law fit Q 2 \F(Q 2 )\ = A(Q 2 /10 GeV 2 )^ 
to the BABAR data with A = 0.182, /3 = 0.25 0. (ii) 
The lower (blue) solid line is also a power-law fit to 
the Belle data with A = 0.169 GeV, fi = 0.18 0], 
whereas the dashed (blue) line below denotes a dipole 
fit Q 2 \F{Q 2 )\ = BQ 2 /(Q 2 + C)to the Belle data with 
B = 0.209 GeV, C = 2.2 GeV 2 Q. (hi) The thick dashed 
line at the bottom shows the dipole fit to the CLEO data 
13 with B = 0.171 GeV, C = 0.6 GeV 2 . 



48| ) that one should not dis- 



tinguish the data as we did in Fig. [5] because the rela- 



The appearance of the Belle data [l[ on the tt° TFF 
forces us to consider the two-photon processes of light 
pscudoscalar mesons in a greater perspective because the 
trend of these data shows an antithetic behavior relative 
to that reported by BABAR 0] with a relative deviation 
of about 1.5a [l[. Instead of a pronounced rise with Q 2 , 
it levels off and follows more or less the scaling behavior 
predicted by QCD and collinear factorization. Though 
this is welcome from the theoretical point of view, an in- 
creasing TFF behavior with Q 2 cannot be ruled out at 
present — at least as long as no possible sources of errors 
have been identified by the BABAR Collaboration to re- 
voke the validity of their results. However, given that 
no unique QCD mechanism has been proposed to pro- 
vide the necessary enhancement of the ir° TFF in order 
to reconcile the BABAR data with QCD, it seems reason- 
able to consider the auxetic behavior of the BABAR data 



5 



and the entailed discrepancy to QCD as an anomaly. 
This view is strengthened by the fact that the pivotal 
QCD effects, notably, the NLO and (the main) NNLO ra- 
diative corrections, ERBL evolution, and the Twist-four 
term give suppression to the TFF, except the Twist-six 
correction which is either small or of the same size as 
the NNLO term but with opposite sign, hence canceling 
against it. Moreover, including in the theoretical analysis 
a small but finite virtuality of the quasi-real photon, as in 
real single -tag ged experiments, also yields to suppression 

0, m m, \m > albeit this 

suppression is more significant 

at lower Q 2 values. 

Bottom line is that skepticism about the accuracy of 
the BABAR data at high Q 2 prevails, giving preference to 
the Belle results that found no deviation from the stan- 
dard QCD scheme. Moreover, in that case there is no 
chasm between the TFFs of the ir° and the \n), i.e., no 
evidence for a significant flavor-symmetry breaking in the 
pscudoscalar meson sector of QCD — in accordance with 
all previous experiments. 
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